Bubble Relaxation Dynamics in Double-Stranded DNA 
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discrete and finite size model for the description of the 
dynamics of bubbles while they used a continuous and 
semi-infinite version of the model to fit their experimen- 
tal data. In this continuous and semi-infinite limit, the 
survival probability of the bubble reads [TJ : 
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The paper deals with the two-state (opening-closing of base pairs) model used to describe the 
fluctuation dynamics of a single bubble formation. We present an exact solution for the discrete 
and finite size version of the model that includes end effects and derive analytic expressions of the 
correlation function, survival probability and lifetimes for the bubble relaxation dynamics. It is 
shown that the continuous and semi-infinite limit of the model becomes a good approximation to 
exact result when a^ <^ 1, where A^ is bubble size and a, the ratio of opening to closing rates of 
base pairs, is the control parameter of DNA melting. 
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Upon heating, a double stranded DNA (ds DNA) un- 
dergoes a denaturation process with the formation of 
bubbles of increasing size and number and, eventually, 
leading to the separation of the two strands [l| . On the 
other hand, many of DNA biological activities require the 
unzipping of the two strands by breaking hydrogen bonds 
between base pairs. Such open regions of complex DNA, 
enclosing up to 10 — 30 broken base pairs, represent a 
first step of the transcription processes and are called the 
transcription bubbles. Several theoretical models have 
been proposed to describe the phenomenon of bubble for- 
mation (for a review see e.g., Q)- However, the issue 
remains unsettled with various, and even contradictory, 
results reported in the literature. This is indicative of the 
complexity of the problem which involves number of fac- 
tors (e.g., base pair sequences, molecular environment, 
counterions, and so on) that can influence the denatu- 
ration process in various ways (see e.g., 0i S 0)- ^^i 
addition, as an one or quasi-one dimensional system, the 
ds DNA is expected to be very sensitive to thermal fluc- 
tuations. Therefore, it seems appropriate in a first step 
to study the fluctuations of local breathing or unzipping 
of a ds DNA that opens up bubbles of a few tens of base 
pairs. 

The characteristic dynamics of these local denatura- 
tion zones (bubbles) in the structure of a ds DNA have 
been recently probed through fluorescence correlation 
spectroscopy [a, Q ■ This is an essential issue not only 
for physiological processes involving ds DNA but also for 
providing insights on the general nature of fluctuations in 
such systems. From a theoretical modeling perspective, 
however, we have just begun to understand these exper- 
imental results. In their recent paper |3, Altan-Bonnet, 
Libchaber and Krichevsky (ALK) have presented a mea- 
surement of the dynamics of a single bubble formation 
in ds DNA construct. The authors proposed a simple 
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where x = t/Tao,c and the bubble lifetime is, 
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where k^ and k- are the opening and closing rates of 
base-pair, respectively, e the bubble extension energy and 
keT the thermal energy. In the same spirit, the dynamics 
of bubble formation have been studied in terms of Fokker- 
Planck equation j^. In this paper, we go one step forward 
in providing the exact solution of the generalized ALK 
model, taking into account both the discreteness of the 
system and the finite size and including end effects. Our 
motivation in this investigation is to provide analytic ex- 
pressions for bubble relaxation function, relaxation time, 
and lifetime. Such exact solutions may significantly im- 
prove data analyzes and be very relevant for any systems 
with arbitrary e and size N . 

Following ALK, we denote by 6,1 (i) the probability 
density of bubbles of size n at time t in the system. As- 
suming that all conformations of the ds DNA can be de- 
scribed as two states (closed or open), the fluctuations 
dynamics in the number n of open base-pairs in the bub- 



ble is described by the master equation, 
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where, in addition to the rates k± in ALK model ll|, we 
have expUcitly introduced the opening and the closing 
rates fci and fc2, respectively, for opening the first and 
closing the last pairs since two ends of the DNA helix are 
sealed. 

Stationary Distribution: When fci 7^ and fc2 7^ 0, 
Eq.Q admits a stationary solution given by. 
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where bst{0) = 1/ [l + {kiQ/k^) + (/cia^/^s)] with Q = 
(1 — a^)/(l — a). The equilibrium fraction of DNA 
molecules that are closed, open and with bubbles in the 
system are given by 6st(0), bgtiN), and /{,, respectively, 
where 
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The equilibrium constants Ki and K2 for the concentra- 
tions of species in the reactions in Fig. ?? are: 



Ki 



[bubble] fci 



and K2 



[open] 
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[closed] fcb [bubble] 

where the backward fcb and forward fc/ rates are. 
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When fci = fc2 = 0, the concentration of bubbles tends 
zero and we have [open] /[closed] = a^ . 

Relaxation Function: To study the fluctuations 
of bubbles, we consider n(z,t|no) = X]n=o ^^ bn{t\no) 
(where &„(i|no) is conditional the probability density of 
finding a DNA molecule with a bubble of size n at time 
t given that the size was uq at time t — 0) the character- 
istic function for the system prepared with the initial 



condition, 6„(t = 0|no) 



The Laplace trans- 



form [n(z,s|no) — /n dt Il{z , t\no) e ***] of n(2,i|no) is 
obtained as, 
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where D(z,s) — fc_|_[z — zi(s)][z — 2:2(5)] and 21^2(5) — 
s/k^ + 1 + a =F y/{s/k- + 1 + a)^ - 4a /2a. The func- 
tions 6o(s|no) and 6jv+i(s|no), obtained by requiring 
that the numerator of n(z, s\no) cancels at the roots of 
D{z, s), are given by, 
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and 
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with 
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fci(l-0i) [szi - k2{l - Zi)] z^ 
fcl(l-Z2) [sZ2~k2{l- Z2)]z^ 
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To fit with the experimental conditions by ALK, we 
assume that the system is prepared in the initial condi- 
tions bst{no)/fb for 1 < no < iV and zero otherwise. The 
quantity of interest is the correlation function CN^t) that 
describes fluctuations in the bubble population at equi- 
librium and is measured by fluorescence correlation spec- 
troscopy method |7|, 
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in which 6„(0|7io) — <^n,no, and we have used the conser- 
vation of the probability density, X]n=o b,i{t\no) = 1. 
Note that Cn{0) = I since 6„(oo|no) ~ bst{n) and 
Cat (00) = 0. Performing the summation in Ea. H12l) . we 
find the Laplace transform of CN{t) as, 
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where F{z) = 

(1 - Z^)[s(l-Z^+l) + (l-^)(fci+fc2Z^)]/[^^A]. 

From this, the bubble relaxation time is obtained as 
Tjv = C{s — 0). Two limiting cases are considered 
depending on fci and fc2. 

• (fci+fc2) > limit: In this case, the bubble relaxation 
time is given by. 
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where H is Heaviside step function defined as H(z) = 
for z < and H(z) = 1 for z > 0. When either fci or 
k2 tends to zero, rjv hnearly decreases respectively with 
either fci or fc2 towards tn(0) defined as, 
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Note that tn{0) is independent of fci and fc2 because 
the kinetics in these hmits is dominated by the bubbles 
decay. As A^ -^ oo, the fluctuations of bubbles become 
independent of N with the relaxation function, 



^ 1 _ fc_(l-a)(l^zi) 
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and lifetime, 
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• fci = fc2 = limit: In this case, 0^(3) = Bn(s), 
where Bpf{t) is the survival probability of bubbles. Like- 
wise, the bubble lifetime r^r — Bn{s — 0) is given by. 



rjv = 



(1 - a^)(l - a^+^) - N{N + 2)(1 - a)^a^ 
fc_(l-a)2(l-a^)(l-a^+i) 



.(18) 



When N ^ 00, Eq.lO reduces to -§00(5) = (1/s) - 
fc_(l — a) [1 — 2i]/s^, and. 
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where Ii[- • •] is the modified Bessel function of order one, 
y = 2xy/a/{l — a)^ and x = t/Toc- It is worth noting that 
even in the N ^ 00 limit the exact solution Ea. (|19|l for 
the bubble survival probability is different from Eq.lQ]) 
given in [3. The fact is that, depending on the size N 
and the parameter "a", the discreteness of the system 
is an ingredient which might be taken into account to 
capture the correct bubble dynamics. This is illustrated 
in Fig. n where the exact survival probability is compared 
with its iV ^ CX3 limit and the ALK continuous model. 
Figure|21shows the departure in the bubble lifetime to the 
continuous limit as a function of bubble size. It clearly 
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FIG. 1: Bubble survival probabilities, from the top to the bot- 
tom, -B]v(t) (solid line), Boc{t) (long-dashed line) and Bao,c{t) 
(dotted line) versus the rescaled lag times I/tm, tJTao and 
t/Toa,c, respectively. 
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FIG. 2: Reduced lifetime, tn/to^ in Eg.llTIll for fci = 
(dashed line) and Eg. 11811 (solid line), as a function of bub- 
ble size, A'^. Quoted numbers represent the bubble extension 
energy e/ksT. 



appears from Figs. ^ and [3 that the continuous limit as 
done by ALK [^ becomes a fairly good approximation 
to exact result for a^ ^ 1 (where a < 1 is the control 
parameter for the ds DNA melting |3, |9(). 

Simple inspection of expressions in Eas. (|13|l . (|14|) and 
((TH|l . and of the figures, indicate that the behavior of 
bubble dynamics is controlled by the bubble size A^ and 
the parameter a (ratio of opening to closing rates of base- 
pairs). As a < 1 according to the experimental situation 
in 7] , the closing of bubbles is the fastest process in the 
bubbles kinetics. The parameter a also controls the de- 
naturation transition. As a — > 1, there is a kind of "crit- 
ical slowing down" where the fluctuations of bubbles are 
described by an unbiased diffusion process. For instance, 
the bubble lifetime in Ea. H18|) reduces to. 
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in the a ^> 1 limit, and t^ diverges with the bubble size. 



TABLE I: Estimate of fc_ using the expressions of the bubble 
lifetime in the case of ki — k2 = 0. In Ref. 7], the experimen- 
tal bubble lifetime is equal to 95 /is at T=303K for TV = 18 
and DNA samples Mig and Ais- 
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It may be useful for practical purposes to have an idea 
of numerical values of physical parameters entering in 
the problem. In the absence of direct measurement of 
k-, for instance, one can use the experimental data in 
[3 in conjunction with theoretical results to estimate the 
closing rate k-. The results of such an estimation are 
presented in Table Q] 

To summarize, we have presented an exact solution 
of the discrete and finite size model in Eq.® for the 
description of the fluctuations dynamics of bubble for- 
mation. The twofold merit of this two-state (open and 
closed) model is to already include sufficient complexity 
of the bubble dynamics over biomolecular relevant scales 
and to allow exact analytical solution. The mains results 
of the paper are the expressions in Eas. H13|l . (|14|l and l|18|) 
for the bubble correlation function, relaxation time and 
bubble lifetime, respectively. These results, consistent 
with available data, may prove to be useful for analysis 
and interpretation of experimental data on bubble fluc- 
tuations and they are amenable for further experimental 
tests. It is worthwhile to mention in addition that dif- 
ferent expressions for the relaxation function and time 
can be generated within the theoretical framework de- 
veloped above by simply using different initial conditions 
in Eq. H12() for the preparation of the system. 

Given the closing and opening rates of base pair, the 
model discussed above allows also to study phenomena 



related to the denaturation mechanisms of DNA such as 
heating, changing buffer surrounding, or applying exter- 
nal torques or forces [ifll, lul llJ, UJI • Likewise, the model 
can easily modified to include more than two states in 
order to describe, for instance, the intermediates states 
between bond and broken states. Finally, although the 
calculations may become more involve and intricate, the 
theory outlined above can be extended in several direc- 
tions in including in Eq.J^J, for example, the effects of 
base pair sequence in the opening and closing rates (two 
and three hydrogen bonds being involved in A-T and G- 
C base pairs, respectively), initiation of several bubbles, 
bubbles fission and fusion processes, and so on. 
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